Mon. Not. R. Astron. Soc. 000, ??-?? () Printed 10 March 2011 (MN I^TeX style file v2.2) 



Recovering cores and cusps in dark matter haloes using 
mock velocity field observations 



■ Rachel Kuzio de Naray^"*"! and Tobias Kaufmann^J 

^ Center for Cosmology, Department of Physics and Astronomy, University of California, Irvine, CA 92697-4575 
. ^ Institute of Astronomy, ETH Zurich- Honggerberg, CH-8093 Zurich, Switzerland 



10 March 2011 



ABSTRACT 

We present mock DensePak Integral Field Unit (IFU) velocity fields, rotation curves, 
and halo fits for disc galaxies formed in spherical and triaxial cuspy dark matter haloes, 
and spherical cored dark matter haloes. The simulated galaxies are "observed" under 
a variety of realistic conditions to determine how well the underlying dark matter halo 
can be recovered and to test the hypothesis that cuspy haloes can be mistaken for 
cored haloes. We find that the appearance of the velocity field is distinctly different 
depending on the underlying halo type. We also find that we can successfully recover 
the parameters of the underlying dark matter halo. Cuspy haloes appear cuspy in the 
data and cored haloes appear cored. Our results suggest that the cores observed using 
high-resolution velocity fields in real dark matter-dominated galaxies are genuine and 
cannot be ascribed to systematic errors, halo triaxiality, or non-circular motions. 

Key words: dark matter — galaxies: kinematics and dynamics — formation — 
hydrodynamics — methods: numerical — methods: N-body simulations. 



. 1 INTRODUCTION 

, Cold Dark Matter (CDM) is one of the main constituents 
' of the current standard model in cosmology, ACDM. N- 
body simulations based on the CDM paradigm predict dark 
. matter haloes best-described by a steep ( "cu spy", "NEW") 
' power law mass density distribution (e.g. iNavarro et al.l 
] Il996al l . While the exact slope of the inner density dis- 
I tribution remains under debate (e.g. iMoore et all 'l999|; 
ICraham etHll 20061 : Istadel et al.llioollNavarr'oe t al. 2010), 
simulations all agree on producing a cuspy profile with a 
slope of a ^ ~1 central part of the halo (see the 

review of lde BlokllioTol ). 

Observed rotation curves of galaxies typically need 
a massive dark matter halo with a nearly constant den- 
sity core in order to describe the data well. Ideal ob- 
jects for observationally studying the "cusp-core problem" 
are low surface brightness (LSB) galaxies which are late- 
type, gas-rich, dark matter-dominated di sc galaxies. Obser- 
vatio ns using long-slit spectra in Ha (e.g. Ide Blok Bosmal 
I2OO2I) or high-resolution two-dim ensional velocity fields 
(|Kuzio de Narav et al.ll2006l . I2OO8I . hereafter K06 and K08, 
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respectively) show that a cored dark matter distribution (de- 
scribed e.g. by an pseudoisothermal sphere profile) provide 
a better fit to the data than a cuspy NEW dark matter pro- 
file. In the few cases where NEW models do reasonably fit 
the rotation curve data, the predicted values for concentra- 
tions C200 are very low and the circular velocity V200 at the 
virial radius very high, inconsistent wit h the cosmological 
(c2oo, V200) — relation, (see the review of lde B lok 2010). 

Most observational problems such as pointing errors, 
centering offsets and non-circular motions are likely too 
small to cause cusps to be mistaken for cores in the observa- 
tions (ide Blokll201Ql l. but halo triaxiality has been claimed 
to explain the prese nce o f cores. The simulation works of 
iHavashi et all (l2007l ) and iBailin et al.1 (|2007l l. for example, 
showed that due to systematic non-circular motions intro- 
duced by a triaxal halo, cusps could be interpreted for cores 
if long-slit rotation curves were used. Similarly, bars induc- 
ing non-circular motions combined with projection effects 
can create the illusio n of a const ant- density core in a circu- 
lar velocity analysis (jValenzu ela et aD l2007l V Pressure sup- 
port in the rotating gas of low mass galaxies could mask the 
underlying dark matter distri bution, but the importan ce of 
this effect remains uncertain (jDalcanton StilDll2010l V 

The aim of this work is to investigate whether cusps in 
(LSB) dark matter haloes could still be mistaken for cores 
if high-resolution t wo-dimensional velocity fie lds data are 
available (as e.g. in lKuzio de Narav et al.ll2009l V In order to 
do so, we use self- consistent high-resolution simulations of 
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(LSB) galaxy formation and "observe" the simulated galax- 
ies using the DensePak Integral Field Unit. We compare the 
mock velocity fields for different signatures from cuspy and 
cored haloes, spherical and triaxial potentials, and the effect 
of supernovae feedback. We also determine how well the un- 
derlying dark matter halo potential can be recovered. 

The outline of the paper is as follows. In Section 2, 
we present the modeling of our initial conditions, the code 
used for the time evolution and the halo parameters derived 
from the simulations. In Section 3, we discuss the process 
of observing the simulations. The mock velocity fields and 
rotation curves are presented in Section 4. In Section 5, we 
present the halo fits to the mock data. A discussion and 
summary are presented in Section 6. 



2 SMOOTHED PARTICLE HYDRODYNAMIC 
SIMULATIONS 

We perform non-cosmological high-resolution N-Body/SPH 
(smoothed particle hydrodynamics) simulations of disc 
galaxy formation by cooling a rotating gaseous mass dis- 
tribution inside equilibrium cuspy spherical, cored spheri- 
cal, and cus py triaxial dark matter haloes similar to what is 
described in lKaufmann et al. (2QQ 7a). The halo parameters 
were chosen to form extendend discs with masses of a few 
X 10^ Mq without bars or massive bulges in a d ark h alo with 
C200 ^ 8-5, as is common in low mass galaxies (|K06[ ). While 
these simulations are not modeling the full complexity of 
disc formation (e.g. mergers are missing) the dynamics of a 
disc dominated by a dark halo (i.e. the late phase of disc for- 
mation) is followed self-consistently and with high numerical 
resolution. 



2.1 Initial conditions 

Our initial conditions for the spherical runs comprise an 
isolated equilibrium cuspy or cored galaxy-sized halo with 
an embedded spinning gaseous mass distribution initially 
in hydrostatic equilibriurn . The halo models are built as 
in iKazantzidis et all (|2004h , and hence they include a self- 
consistent description of the velocity distribution function. 
Those models with virial mass of M200 = 1.3 x 10^^ Mq 
were constructed using the a/37 models as described in 
IZemp et al.l (|2008l ) . Cuspy haloes were modeled using the 
NFW form (1, 3, 1 for a,/3,7) and C200 = 8.5 following 
iMaccio et al.l (|2008l ). For cored haloes, the parameters were 
set to 0.8, 3, and C200 = 24. This choice produced an ini- 
tial core with a radius of ^ 2 kpc. similar to what has been 
found in observational work (|K06r ). 

We initialise a fraction of the virial halo mass, /b = 0.01, 
as a hot baryonic component with the same radial distribu- 
tion as the dark matter. The initial temperature structure 
of the gaseous halo is calculated by solving the equation for 
the hydrostatic balance of an ideal gas inside a dark mat- 
ter halo. The gas halo is spun with a gas spin parameter of 
Xg — 0.06, and the dark matter particles are init i alised w ith 
no net angular momentum fsee lKaufmann et al.l (|2007al ) for 
details). We note that the baryon fraction is much lower 
than the universal baryon fraction o f ^ 0.17, though this is 
not unexpected in low mass galaxies (jMcGaugh et al.ir2010h . 
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Figure 1. The ratio between the short and long axes of the 
dark matter density distribution of the triaxial halo before cool- 
ing (black line) and at the time of the disc analysis (red line). 



This choice of /b forms discs with masses of a few xlO^M© 
in our models. 

The models use 2.5 x 10^ dark matter and 2 x 10^ gas 
particles (mass of a dark particle ^ 7 x 10^ Mq) and have a 
spatial resolution (softening length) of 100 pc. 

To construct an equilibrium triaxial halo in which the 
gas acquires its angjular momentum throug;h an equal mass 
merger fsee iMoore et al.ll20o3 : iKaufmann et al]l2007al ). we 
start with two spherical, non-rotating NFW haloes each with 
half of the mass and particles of the fiducial model consid- 
ered previously. The two haloes are placed at a separation 
of twice their virial radius. One of them is given a transverse 
velocity of 27 km s~^, and the haloes are allowed to merge. 
This net velocity (determined using a trial-and-error proce- 
dure) results in a final gas distribution that has a similar 
spin parameter to that used in the spherical models. The 
triaxial halo resulting after the merger is first evolved for 
15 Gyr with an adiabatic equation of state for the gas. The 
inner part of the dark halo is significantly triaxial, c/a ~ 0.6 
(and b/a ^ 0.7) and the subsequent formation of the disc 
due to cooling modifies the s hape of the dark matter only 
weakly (see Figure [T] and also IKazantzidis et al.l (|2010[ )V 

2.2 Time evolution and derived halo parameters 

All of the systems were then evolve d for ^ 4.5 G yr using the 
parallel TreeSPH code Gasoline (|Wadslevll20Q4l ). which is 
an extension of the pure N -Body gravity code Pkdgrav 
developed by IStadell (|200lh . The code implements Comp- 
ton and radiative cooling, as well as the star for mation and 
''blast wave feedback" pre s cripti ons detailed in iKatz et al.l 
(|l996h and IStinson et al.1 (|2006l ). In order for star parti- 
cles to form, the gas temperature must be < 30, 000 K, 
the local gas density must be > 0.1 cm~^ and other lo- 
cal gas particles must define a converging flow. Stars then 
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Table 1. Dark halo parameters derived from the simulations after 
the evolution with cooling 



Name 


R2OO or Tcore 


C2OO or pcore 


Cuspy 


104.8 kpc 


8.4 


Cuspy SN 


104.3 kpc 


8.3 


Tri axial 


102.9 kpc 


8.8 


Cored 


1.5 kpc 


69 X lO-^M© pc-3 



form according to th e Schmidt law (j Kennicuttll 1 9981 ) using a 
iMiller Scalol (| 19791) initial mass function. The supernova 
(SN) feedback model creates turbulent motions in nearby 
gas particles that keep them from cooling and forming stars. 
Because of the lack of metal and molecular cooling, the ef- 
ficiency of our cooling function drops rapidly below 10^ K 
ac ting as a temperature flo or similar to what has been used 
in lKaufmann et aL l (|2007bl l. Star formation is only switched 
on in the cuspy model using the feedback prescript iorQ; the 
other simulations are cooling-only runs. 

At the time of the analysis, thin, extended baryonic 
discs resembling LSB galaxies have formed in the centre of 
the haloes. This redistribution of the baryonic matter due 
to cooling also affects the underlying dark matter distri- 
bution, see below. The discs do not show bars or massive 
bulges, have baryonic masses of 5 x lO^M©, and star for- 
mation rates, SFR ^ Q.33 Mc;) yr~^, conapara ble to those 
observed in LSB galaxies (Wyder et alJ I2OO9I ). The discs 
are dark matter-dominated with the baryonic mass reaching 
up to 30% of the total mass within the inner 0.5 kpc and 
far less beyond. The galaxies formed in the spherical sim- 
ulations lie between the mean and the upper 2a deviation 
of the (observationa l) Radius- Velocity relation presented in 
iDutton et al.l (|2007l l, while the galaxy formed in the triax- 
ial simulation lies just above the upper 2a deviation from 
the mean. These somewhat large radial scale-lengths are not 
unexpected as the spin parameter chosen in the initial con- 
ditions was above the m ean. A dditionally, it has been sug- 
gested by Mac cio et al] (| 20071 ) that LSB galaxies reside in 
haloes with spins higher than average. 

We derive halo parameters after the evolution with cool- 
ing directly from the simulations and list them in Table (TJ 
Additionally, in Table [1] we establish our naming conven- 
tion for the simulations. We refer to the spherical cuspy 
simulation without feedback as 'Cuspy'; the spherical cuspy 
simulation with feedback is referred to as 'Cuspy SN'. The 
triaxial cuspy simulation is referred to as 'Triaxial' and the 
spherical cored simulation is referred to as 'Cored'. Moti- 
vate d by the param eterisations used in obser vational works 
fe.g. Ide Blokll2010l ) we fit an NFW profile (|Navarro et al.1 
Il996af r to the dark matter density curves derived from the 
cuspy simulations and a pseudoisothermal sphere profile (as 
in e.g. |k06|) to the cored simulation. 

The initial and final halo shapes do not differ signifi- 
cantly, even in the simulation with feedback. In our simula- 
tions, massive dark matter haloes are in place at the begin- 
ning of the simulation and the low mass present in baryons 



^ Heating by a unifor m UV background from QSO following 
iHaardt Madaul {1996) has also been used. 



{fbaryons = 0.01) is of little importance for the dynamics of 
the dark matter (specifically, the mass of blown-out baryons 
is too small to affect the dynamics of the dark halo). The 
only potential infiuence from the supernova feedback is a 
disturbance of the velocity field of the gas in the galactic 
disc. This is diff e rent t han the cosmological simulations of 
iGovernato et al ] (|2010l ) (using the same feedback prescrip- 
tion but a higher minimal density for star formation) who 
find that cored dark matter haloes are created due to strong 
outfiows from super novae in dwarf galaxies. The outfiows 
remove the central dark matter cusp at high redshift when 
the dark matter halo is less massive. 



3 "OBSERVING" THE SIMULATIONS 

Each galaxy+halo simulation is "observed" under a va- 
riety of instrumental configurations and physical proper- 
ties designed to represent a range of mock data quality 
from ideal to realistic. We test three distances, 15 Mpc, 
30 Mpc, and 45 Mpc, corresponding to spatial resolutions 
of - 0.07 kpc/'', 0.15 kpc/'', and 0.22 kpc/", respectively. 
The choice of these spatial resolutions was motivated by the 
spatial resolution obtained for the LSB galaxies observed 
in K06 and IKOSI , as well as a requirement that a 3" IFU 
fiber provide sub-kiloparsec resolution. We test three galaxy 
inclinations, i = 30°, 50°, 65°, th e ran ge again being mo- 
tivated by the data presented in IK06| and IKOsI . There is 
an inherent trade-off between inclination and obtaining a 
well-sampled velocity field, and most observational stud- 
ies target galaxies having inclinations within, or very close 
to, th is range (e.g. ISpano et al.ll2008l : iBlais-Ouellette et al.l 
12004 ). Not all lines of sight in nonaxisym metric systems are 
equivalent ('K uzio de Narav et al ] 12009'). so we also rotate 
the galaxy+halo system by a random angle in the plane of 
the disc. For consistency, we do this for the spherical simu- 
lations in addition to the triaxial simulation. 

DensePak, an integral field spectrograph on the 3.5- 
m WIYN telescope at the Kitt Peak National Observa- 
tory (KPNO), is used to "observe" each simulated galaxy. 
DensePak is a 43" x 28" fixed array of 3" fibers with 3.84" 
separations. We model the 85 working, and five missing or 
broken, fibers in the main bundle. We test three configu- 
rations of DensePak fiber alignment (see Figure [2]) . In the 
Ideal case, 11 pointings of the IFU array provide maximum 
spatial coverage across the galaxy. In the Realistic 1 and 
Realistic 2 cases, there are 5 pointings of the IFU array. 
In the Realistic 1 case, the IFU fibers cover the central re- 
gion of the mock galaxy. The IFU fibers are along the mock 
galaxy major axis in the Realistic 2 configuration. The ve- 
locity recorded for each fiber is the average of all the line-of- 
sight velocities of the gas particles falling within the fiber. 
In the Ideal case, all 935 IFU fibers provide a velocity mea- 
surement. To simulate the sparse Ha emission observed in 
some real galaxies, ^ 40% of the fibers in the Realistic 1 and 
2 cases are randomly turned off so that no velocity measure- 
ment is obtained. Finally, we ensure that the fiber-to-fiber 
velocity variation is 6-10 k m s~^ , consistent with the velocity 
dispersion measured in the lKOd galaxies. 

Mock velocity fields are obtained using all three IFU 
fiber configurations at the high and medium spatial resolu- 
tions. Only the Ideal and Realistic 2 configurations are used 
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at low spatial resolution because of the reduced angular size 
of the galaxy. The three galaxy inclinations are tested at 
each spatial resolution and IFU fiber configuration. For each 
galaxy+halo simulation, we obtain five mock velocity fields 
at each combination of galaxy parameters and IFU configu- 
rations. Our final data set includes 120 mock velocity fields 
per galaxy+halo simulation, for a total of 480 mock velocity 
fields. 



4 MOCK DATA 

In this section, we present examples of the mock veloc- 
ity fields and derived rotation curves for each of the four 
galaxy+halo simulations. 
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Figure 2. Configurations of the IFU fibers (red circles) shown on 
the galaxy (grey points). Top: There is maximum spatial coverage 
across the galaxy in the Ideal configuration and all fibers provide 
a velocity measurement. Middle: In the Realistic 1 configuration, 
the reduced number of fibers cover the central part of the galaxy. 
Additionally, ~40% of the fibers shown are randomly turned off so 
that not all fibers provide a velocity measurement. Bottom: The 
fibers in the Realistic 2 configuration are along the galaxy major 
axis. As in the Realistic 1 configuration, not all fibers provide a 
velocity measurement. 



4.1 Mock velocity fields 

In Figures [3]|6l we show four representative examples of the 
mock velocity fields obtained for each simulation, covering 
the range of spatial resolutions, inclinations, and IFU config- 
urations. The velocity fields in the upper left and lower right 
are examples of the data obtained using the Ideal configura- 
tion while the upper right shows an example of the Realistic 
1 configuration and the lower left an example of the Realis- 
tic 2 configuration. Figure [3] contains mock velocity fields for 
the triaxial simulation. Mock velocity fields for the spherical 
cuspy simulation are in Figure H] those including feedback 
are in Figure [5] Figure [6] contains mock velocity fields for 
the cored simulation. 

In general, the mock velocity fields of the triaxial simu- 
lation (Figure [3]) are characterised by misaligned kinematic 
and photometric major axes, and non-perpendicular kine- 
matic major and minor axes (see Figure [7|). Most notice- 
ably, the minor axis is often twisted. Though no physical 
bar is present in the simulation, these are features com- 
monly obs^vedjn^^ of barred galaxies (e.g. 
lBosmalll98ll : [Hernandez et al.ll2005l ) . Bar-like dynamical sig- 
natures (without an obvious bar) are also expected based 
on simulation work for subdominant discs in triaxial halos 
iKazantzidis et al.ll2010i ) . The orientation of the twist varies 
in the different mock velocity fields depending on the angle 
between the observer's line-of-sight and the potential. 

The kinematic and photometric major axes are aligned 
in the mock velocity fields of the spherical cuspy simulations 
(Figures m and [5|), and the major and minor kinematic axes 
are perpendicular. The isovelocity contours are noticeably 
pinched at the centre of the velocity field. There is very little 
difference between the mock velocity fields of the simulations 
with feedback and those without. The deposition of energy 
due to the supernovae did not affect the average velocity field 
in the cold gas heavily, thus it is still tracing the underlying 
(dark matter) potential quite closely. 

The mock velocity fields of the cored simulation (Fig- 
ure [6]) are characterised by parallel isovelocity contours at 
the centre. The kinematic major and minor axes are perpen- 
dicular, and the kinematic and photometric major axes are 
aligned. The central regions of these velocity fields are very 
similar to those of dwarf and LSB galaxies (e.g. Gentile et 
al. 2005; K06). 

We find through visual inspection of the 480 mock ve- 
locity fields, that depending on the underlying halo type, the 
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Figure 3. Sample mock velocity fields of the triaxial cuspy sim- 
ulations; isovelocity contours are at 10 km s"-*^ intervals. The 
velocity fields in the upper left and lower right are examples of 
the Ideal configuration; the upper right and lower left are exam- 
ples of the Realistic 1 and Realistic 2 configurations, respectively. 
The kinematic and photometric major axes are not aligned, the 
kinematic major and minor axes are not perpendicular, and the 
minor axis is often twisted. 




Figure 5. Sample mock velocity fields of the spherical cuspy SN 
simulations; isovelocity contours are at 10 km s"-*^ intervals. The 
layout is the same as Figure [3] The kinematic and photometric 
major axes are aligned, the kinematic major and minor axes are 
perpendicular, and the isovelocity contours are pinched at the 
centre of the velocity field. 




Figure 4. Sample mock velocity fields of the spherical cuspy 
simulations; isovelocity contours are at 10 km s"-*^ intervals. The 
layout is the same as Figure [3] The kinematic and photometric 
major axes are aligned, the kinematic major and minor axes are 
perpendicular, and the isovelocity contours are pinched at the 
centre of the velocity field. 




[kmf<} 

Figure 6. Sample mock velocity fields of the spherical cored simu- 
lations; isovelocity contours are at 10 km s"-*^ intervals. The layout 
is the same as Figure [3] The kinematic and photometric major 
axes are aligned, the kinematic major and minor axes are perpen- 
dicular, and the isovelocity contours are parallel at the centre of 
the velocity field. 



mock velocity fields look distinctly different. We capture the 
distinguishing features discussed above in the illustration in 
Figure El We can quantify these differences by comparing 
velocities measured along slits placed parallel to the minor 
axis. In Figure [9l we plot the average of the measured ve- 
locities in a 3^^ x 28" slit (the width of a DensePak fiber 
and the width of the DensePak array, respectively) parallel 
to and offset from the minor axis by ±0.5 kpc. The velocity 
is roughly constant as a function of position along the slit 
in the cored mock velocity fields, consistent with the gener- 
ally parallel isovelocity contours. The pinch observed in the 
spherical cuspy mock velocity fields translates to a symmet- 
ric ^ 8 km s~^ rise in velocity along the slit as the major 



axis (r = 0) is approached. When slits are placed parallel 
to the photometric minor axis in the triaxial mock velocity 
fields, there is an asymmetric (around r = 0) overall change 
in velocity of ^ 16 km s~^ along the slit. This behavior is 
unique to the triaxial case; velocities measured along slits 
placed at an angle to the minor axis in the spherical cuspy 
and cored cases show a steady change (constant increase or 
decrease) in velocity along the entire length of the slit. It is 
important to note that the differences in the central regions 
of the mock velocity fields are equally detectable at all three 
spatial resolutions. That the underlying halo type produces 
a unique and measurable signature in the observed velocity 
field is an important result. It means that the halo profile 
can be probed in a model-independent way. 
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Figure 7. Difference in degrees between the photometric and 
kinematic major axes as measured at a radius of 19.5^^, approxi- 
mately half the length of a single DensePak pointing. The axes are 
aligned in the majority of the spherical cuspy (with and without 
feedback) and cored mock velocity fields; the axes are frequently 
misaligned in the triaxial mock velocity fields. 
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Figure 8. Schematic representations of the distinguishing fea- 
tures of the triaxial cuspy, spherical cuspy, and spherical cored 
mock velocity fields. 

4.2 Application to LSB galaxy data 

We apply the slit method described in Section 14.11 to the 
Dens ePa k vel ocity fields of the LSB galaxies presented in 
iKOd and IKOsI to see how well it compares to the results ob- 
tained for the galaxies by fitting halo models to the rotation 
cur ves. W e p lot th e results for 5 of the galaxies in Figure 
[TOl IKOGI and IKOSI find the cuspy NFW halo to provide a 



Figure 9. Average velocity measured along slits parallel to, and 
offset from, the minor axis of the mock velocity fields. A typical 
errorbar is shown in the lower left corner. 

better fit to the rotation curve of NGC 4395 than the cored 
isothermal halo; they also note a misalignment of the minor 
axis. In Figure IIOI the change in velocity along the slit in 
NGC 4395 displays the same shape as those along the triax- 
ial slits in Figure [9l There is an asymmetric (around r = 0) 
overall change in velocity of ^ 18 km s~^ along the slit. 
Within the errors, the velocities are roughly constant along 
the slits for UGC 4325, DDO 64, and F583-1. A constant ve- 
locity is expected for cored haloes, and these three galaxies 
are found t o be best-described by cored isothermal haloes 
in K06 andlKoi. [k08 find the rotation curve of NGC 7137 
to be equally fit by cored and cuspy haloes. In Figure IIOI 
the velocities along the slit in NGC 7137 display a close to 
symmetric (around r = 0) ^ 6 km s~^ change in velocity. 
This behavior is similar to the change in the velocity seen 
along the spherical cuspy slits in Figure [9l We find that the 
slit method can provide a measure of the underlying halo in 
real galaxy data that is in agreement with results based on 
the traditional method of rotation curve fitting. 

4.3 Mock rotation curves 

The NEMO (jTeubenl [l995l l program ROTCUR ("Bege manl 
1 19891 ) is used to derive a rotation curve from each mock ve- 
locity field. ROTCUR treats the observed velocity field as 
an ensemble of tilted rings and fits for the centre, systemic 
velocity, inclination, position angle, and rotation velocity in 
each ring. We allow the position angle to change with radius, 
but we do not make any explicit attempts to model the ob- 
vious non-circular motions when fitting the triaxial velocity 
fields. For an extensive explanation of ROTCUR and its 
application to the Den sePak velocity fields, the reader is re- 
ferred tolKOelandlKoi. We also calculate the rotation curve 
directly from each of the four galaxy+halo simulations as- 
suming — GM/r. 

In Figure [TTl we compare the rotation curves extracted 
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Figure 10. Average velocity measured along slits parallel to, 
and offset from, the minor axis of galaxies from K06 and .K08. . A 
typical errorbar is shown in the lower left corner. 



from each mock velocity field to the rotation curves derived 
from the simulations. Overall, we find the mock rotation 
curves to trace the simulation rotation curves well. The ro- 
tation curves extracted from the triaxial velocity fields show 
a large spread in velocity at all radii, not surprising given 
the dependence on the observer's line-of-sight with respect 
to the orientation of the potential. There is much less varia- 
tion in the mock rotation curves extracted from the spherical 
cuspy and cored velocity fields. In addition, the comparisons 
of the average mock rotation curves to the rotation curves 
derived directly from the simulations are practically indis- 
tinguishable. In all three cases, the average mock rotation 
curve is only ^3 km s~^ below the simulation rotation curve 
at small radii ( r < 1 kpc), is equal to the simulation at in- 
termediate radii, and displays a trend toward slightly higher 
average velocities at large radii. 



5 HALO FITS 

Following the convention of observational studies, we fit each 
mock rotation curve with both the cuspy NFW halo and 
the cored pseudoisothermal halo. We fit the data under the 
assumption of negligible baryonic mass. The NFW rotation 
curve is given by 



ViR) = V2i 



ln(l + cx) — cx/(l + cx) 
x[ln(l + c)-c/(l + c)] ' 



fl) 



with X = R/R200. Here R200 is the radius at which the den- 
sity contrast exceeds 200, V200 is the circular velocity at 
i?2oo, and the concentration parameter c = R2oo/Rs, where 
Rs is the ch aracteristic radius of the halo ([Navarro et al.l 
Il996al . I1997I I. The rotation curve of the pseudoisothermal 
halo is given by 
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Figure 11. Comparison of the recovered mock rotation curves 
(black open circles) to the rotation curve derived from the sim- 
ulation (red line) using = GM/r. The green stars are the 
average mock rotation curve. The bottom panels plot the ratio of 
the simulation rotation curve to the average mock rotation curve. 



8 Kuzio de Naray & Kaufmann 



1 



Rc , ( R 

— — arctan -— 
R \Rc 



(2) 



where po is the central density of the halo and Rc is the 
core radius of the halo. 

In Figure [T2l we show the distribution of the recovered 
halo parameters (c, i?2oo; po, Rc) from the fits to the mock 
data (In the following, i?2oo, po and Rc are given in units 
of kpc, 1O~^M0 pc~^ and kpc, respectively). We find that 
the parameter values obtained from the fits to the mock 
rotation curves are in agreement with those derived directly 
from the simulations (Table [T]). The parameters of the best- 
fitting NFW halo to the spherical cuspy data are c — 9.4itl.8 
and i?2oo = 107 =b 10; for the spherical cuspy with feedback 
data, the values are c = 8.5 di 3.0 and i?2oo = 113 =b 20. 
There is more variation in the halo fits to the triaxial cuspy 
data. The best-fitting halo parameters are c — 11.0 ±4. 5 and 
i?2oo = 87=bl0. The parameters of the best-fitting isothermal 
halo to the spherical cored data are Rc — 2.0 ± 0.5 and 
po = 50 ± 5. 

Fitting a pseudoisothermal sphere to the cuspy data re- 
sults in dense, small cores (po ^80; Rc ^1-0); similarly, when 
fitting cuspy haloes to the cored data, the data want c < 1.0 
and high i?2oo- This is not surprising since a profile with a 
tiny core and high central density resembles a cuspy pro- 
file and only with a very low concentration can one stretch 
an NFW curv e to fit the rotation curve of a cored halo, see 
Ide Blok" (2010'). However, given the information from the ve- 
locity fields, one can easily break this ambiguity and iden- 
tify the underlying halo shape even before fitting a profile 
to it. Additionally, t he mass-concentration relation of CDM 
(jMaccio et al.l [ioosl l would predict c < 1.0 halos to have 
rotation velocities many times higher than the galaxy itself. 

In Figure [13] are representative examples of residual 
mock velocity fields showing the differences between the 
mock data and the best-fitting cuspy and cored halo mod- 
els. For the cored mock velocity fields, we display the resid- 
uals obtained when subtracting the best-fitting spherical 
cuspy halo model; the residuals obtained when comparing 
the cored data to the best-fitting spherical cuspy with feed- 
back halo model are comparable. The residuals are small, 
typically less than ^ 4 km s~^, when cusps are fit to cusps 
and cores are fit to cores. When a core is fit to the cuspy 
data or a cusp is fit to the cored data, the residuals can 
become quite large. Residuals as high as '-^ 15 km s~^ and 
^ 30 km s~^ are obtained when cores are fit to the spherical 
and triaxial cuspy data, respectively; residuals on the order 
of 15 km s~^ are obtained when cusps are fit to the cored 
data. 

We also calculate the inner slope of the mass density 
profile by converting the mock rotation curves using 



p{R) 



1 

4^ 



^Vd_V fV_. 
RdR \r) 



(3) 



where G is the gravitational c onstant and V" is the observed 
rotation velocity at radius R (|de Blok et ahllioOll ). 

In Figure [TH we show the distribution of the asymp- 
totic inner slope of the dark matter density profiles derived 
from the mock rotation curve data. The gray shaded regions 
indicate the inner slopes derived directly from the dark mat- 
ter densities of the simulations. We find the cuspy halos to 
generally have inner slopes of —1 or steeper and the inner 
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Figure 12. Recovered halo parameters from the mock velocity 
fields and rotation curves. The contours represent the distribu- 
tion of parameters found with the mock observations and the 
red crossed circles indicate the values of the parameters derived 
directly from the simulations. 



slope of the cored halo to be more shallow, around a value of 
^ —0.4. The inner slopes derived directly from the simula- 
tions are —1.3, —1.4, —1.4, and —0.3 for the spherical cuspy, 
spherical cuspy with feedback, triaxial cuspy, and spherical 
cored halos, respectively. The corresponding inner slopes de- 
termined from the mock data are — 1.0 ± 0.7, — 1.0 ± 0.8, 
-1.4 ±0.7, and -0.5 ± 0.4. 



6 DISCUSSION AND CONCLUSIONS 

We have presented mock velocity fields, rotation curves, and 
halo fits for simulated LSB galaxies formed in spherical and 
triaxial cuspy dark matter haloes and spherical cored dark 
matter haloes. The mock velocity fields span a range of data 
quality representing ideal to realistic observations. The main 
findings of this work are: 

• The underlying halo type produces a unique signature 
in the velocity field. This can be used to constrain the shape 
of the dark matter profile (spherical or triaxial cuspy, or 
spherical cored) without fitting an analytic density profile 
to it. 

• Cored and cuspy haloes can also be distinguished 
clearly by deriving their asymptotic inner slopes from the 
rotation curve data. 

• Given at least one of the above information one can then 
successfully recover the underlying halo parameters from the 
rotation curve using the appropriate analytic form for the 
density profile (NFW or pseudoisothermal sphere). 

This means that if an LSB galaxy were in a cuspy halo, the 
cusp would he observable in the data. Given these results, 
we find it difficult to mistake cuspy haloes, spherical or tri- 
axial, for cored haloes. The observed cores in dark matter- 
dominated galaxies are true discrepancies from the predic- 
tions of (dark matter-only) ACDM simulations. Systematic 
effects, non-circular motions, and halo triaxiality cannot ex- 
plain the observed differences. 
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Figure 13. Sample residual velocity fields after subtracting the 
best-fitting cuspy and cored halo models from two different ob- 
served mock velocity fields (Realistic 1 case on the left, Realistic 
2 case on the right.) 
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Figure 14. Distribution of the asymptotic inner slope of the dark 
matter density profiles derived from the mock rotation curve data. 
The gray shaded regions indicate the inner slope derived directly 
from the dark matter densities of the simulations. 



Baryonic processes more effective than those we have 
modeled or entirely different dark matter models may be 
necessary to explain the observed cores. Feedback from star 
formation could transform the initial cusps into cores. If 
LSB galaxies form with a massive dark matter halo already 
in place and with the same low baryon fraction that we ob- 
serve today, then feedback from star formation of the type 
we have modeled would not be enough to change a cusp into 
a core. In this scenario, the dark matter halos in which LSB 
galaxies reside cannot be cuspy CDM halos. If, however, 
feedback processes from star formation can affect dark mat- 
ter haloes during the early stages of their formation, it is pos- 
sible that cusps c an be transformed into c ores. Thi s has been 
demonstrated bv lGovernato et al ] (|2Q10l ) (see alsolOh et al.l 
([2010 )) using a set of fully cosmological simulations of dwarf 
galaxy formation {Mvir 3 x 10^° M©) where strong out- 
flows from supernovae remove low-angular-momentum gas 
and decrease the dark matter density to less th an half in the 
central part ( see also Navarro et al. 1996b; Re ad k, Gilmord 
[2OQ5I : tMashche nko et al. 2008). The question in this scenario 
then becomes one of detecting the blown-out baryons. It 
also remains to be seen if this kind of mechanism can act 
efficiently on higher mass scales, given the deeper potential 
wells. 

Alternatively, the observed cores may need to be ad- 
dressed by a different dark matter model/particle that nat- 
urally produces cored haloes. Recent results, however, have 
found that self-interactions and warm dark matter particle 
properties cannot be responsible for the cores observed in 
dark matter-dominated galaxies (e.g. |Kuzio de Narav et al.l 
I2OIOI : IVillaescusa-Navarro Dalalll2010l V 
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